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Coherence effect in a two-band superconductor: Application to iron 

pnictides 

Keisuke Masuda* and Susumu Kurihara 

Department of Physics, Waseda University, Okubo, Shinjuku, Tokyo 169-8555, Japan 

From a theoretical point of view, we propose an experimental method to determine the 

pairing symmetry of iron pnictides. We focus on two kinds of pairing symmetries, s_| and 

s++, which are strong candidates for the pairing symmetry of iron pnictides. For each of 
these two symmetries, we calculate both the density and spin response functions by using 
the two-band BCS model within the one-loop approximation. As a result, a clear difference is 

found between the s_| - and s ++ -wave states in the temperature dependence of the response 

functions at nesting vector Q, which connects the hole and electron Fermi surfaces. We point 
out that this difference comes from the coherence effect in the two-band superconductor. We 
suggest that the pairing symmetry could be clarified by observing the temperature depen- 
dence of both the density and spin structure factors at the nesting vector Q in neutron 
scattering measurements. 

KEYWORDS: iron-pnictide superconductors, coherence effect, pairing symmetry, response 
functions, two-band BCS model 



1. Introduction 

Recently discovered iron pnictides 1 -* are expected to be unconventional (non-phonon- 
mediated) superconductors. For such superconductors, determining the pairing symmetry is 
of great importance, because the pairing symmetry is closely related to the mechanism of 
superconductivity. 

In iron pnictides, many theoretical studies have focused on the pairing symmetry. In the 

early stage, Mazin et al. 2 ^ and Kuroki et al. 3 ^ suggested that the -wave state, in which 

the s-wave gap on the hole Fermi surface has the opposite sign as that on the electron Fermi 
surface, is favored due to the nesting between the hole and electron Fermi surfaces. This 
s-i -wave state has also been supported by other studies; the theory based on the random- 
phase approximation (RPA) using the 16-band d-p model, 4 ) the perturbation theory, 5 ) and 
the functional renormalization-group theory. 6 ' 7 > On the other hand, within the fluctuation 

exchange approximation (FLEX), it has been shown that the s_| -wave state is not necessarily 

favored. 8 ) Onari and Kontani pointed out that the s_| -wave state is quite fragile against 

nonmagnetic impurities, and suggested that the s++-wave state, in which the s-wave gap on 
the hole Fermi surface has the same sign as that on the electron Fermi surface, is a promising 
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Fig. 1. (a) Schematic temperature dependence of the acoustic attenuation rate, (b) Schematic tem- 
perature dependence of the nuclear spin relaxation rate. 



candidate for iron pnictides. 9 ) Recently, by the microscopic calculations, they showed that the 
s ++ -wave state is induced by orbital fluctuations due to the electron-phonon interaction. 10 ) 

Many experiments have also been performed to clarify the character of iron pnictides. The 
disconnected Fermi surfaces (some hole and some electron Fermi surfaces) have been observed 
by angle-resolved photoemission spectroscopy (ARPES). 11,34 ) In addition, the fully gapped 
s-wave nature has been confirmed by ARPES 13 ~ 15 ) and penetration depth measurements. 16 ) 
However, there are few experimental methods which identify the phase difference between 
the superconducting gap on the hole Fermi surface and that on the electron Fermi surface. 
Therefore, it is essential to suggest an experimental method to reveal the relative phase 
between the two superconducting gaps. 

In this paper, we suggest an experimental method which distinguishes between the - 

and s ++ -wave states. In order to find the difference between these two states, we focus on the 
coherence effect in the response functions. The coherence effect is one of the most important 
evidences that determined the success of the BCS theory. 17 ) Experimentally, a qualitative 
difference between the behavior of the acoustic attenuation rate 18 ) and that of the nuclear spin 
relaxation rate 19 -* had been established. The former decreases monotonically as temperature 
decreases [see Fig. 1(a)]. In contrast, the latter has a peak (called the "coherence peak") just 
below T c [see Fig. 1(b)]. Although the phenomenological two-fluid model could not explain 
both behaviors, 20 ^ the BCS theory gave the correct explanation of both features. 

According to the BCS theory, these behaviors can be understood as a result of the interfer- 
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ence effect (called the "coherence effect") of the field-induced fluctuations in superconductors. 
The temperature dependence of the response function is determined by the time reversal sym- 
metry of the perturbing field. 21 > The perturbing field which relates to density fluctuations is 
even under time reversal. As a result, the density response function which corresponds to the 
acoustic attenuation rate has the following coherence factor. 



This coherence factor has a value close to near the Fermi level, and cancels the square- 
root divergence in the density of states (DOS) which occurs in the superconducting state. 
Hence, the density response function shows a monotonic decrease as temperature decreases. 
In contrast, the perturbing field which relates to spin fluctuations is odd under time reversal. 
As a consequence, the spin response function which corresponds to the nuclear spin relaxation 
rate includes the following coherence factor. 



divergence in the DOS. Hence, the spin response function has a peak just below T c . This is 
the theoretical explanation of the coherence effect in the single-band s-wave superconductor. 
It is worth noticing that the opposite behavior of the two different response functions can be 
explained within a framework of the BCS theory. This is the reason why the coherence effect 
has become the strong evidence for the BCS theory. 

To investigate the coherence effect in iron pnictides, let us extend the above analysis to 
multi-band cases. First we introduce the two-band BCS model, and then calculate both the 
density and spin response functions within the one-loop approximation. As mentioned above, 

the coherence factor includes the product of the two gap functions (AkAk+ q ). In the s H - 

wave state, this product has a negative value (A^A^+q < 0) when q equals Q, which is the 
nesting vector connects the hole and electron Fermi surfaces. Then the sign in front of the 
fraction in the coherence factor is effectively reversed. Therefore, in this state, it is expected 
that the density (spin) response function has a peak just below T c (has a monotonic decrease 
as temperature decreases) as opposed to the single-band s-wave case. In fact, these behaviors 
are also confirmed in our calculations. On the other hand, in the s ++ -wave state, the product 
of the two gap functions has a positive value (AkAk+ q > 0). Therefore the sign in front of the 
fraction in the coherence factor does not change, and the behavior of the response functions 
is the same as in the single-band s-wave case. 

From our analysis it follows that the pairing state would be identified by observing the 
temperature dependence of both the density and spin response functions (density and spin 
structure factors in actual experiments) at the nesting vector Q. 



1 / _ AkAk+q \ 

2 \ ^k-Ek+q / 




This coherence factor has a value close to 1 near the Fermi level, and does not cancel the 
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2. Two-band BCS model 

In this section, we explain the two-band BCS model used in the calculation of the response 
functions. As a starting point we use the tight binding model proposed by Raghu et al. 22 ^ 
This model includes two atomic orbitals (d xz and d yz ), and the hopping integrals up to the 
second nearest neighbors. The tight binding Hamiltonian T-Lq is given by 

«o = X>J(k) [Mk)-/*) 1 + e -(k)r 3 + txyWn] Vv(k), (1) 

k<7 

with 

where d xa (k) (d ya (k)) is the annihilation operator with spin a and wave vector k in the orbital 
d X z (dyz)- T~i are the Pauli matrices, and 

€±(k) = i[€ x (k)±€ y (k)], 

e x (k) = — 2t\ cos k x — 2t2 cos k y — 4tz cos k x cos k y , 

e y (k) = — 2*2 cos k x — 2t\ cos fc y — 4tz cos /c x cos k y , 

e xy (k) = — 4*4 sin k x sink y . (3) 

Here, t\ and *2 (£3 and £4) are the nearest-neighbor (next-nearest-neighbor) hopping integrals. 

We introduce the band operator (k) which annihilates an electron with spin a and 
wave vector k in the band v: 

d ra (k) = a r u (k) lua (k), (4) 

v=± 

with 

<(k) = a y _(k) = sgn[e xy (k)] 



1 + c-(k) 



^ 2 2^62 (k) + e y k) - 



<(k) = -a*(k)= 1- € - (k) = . (5) 

\J 2 2^ 2 (k) + e 2 j/(k ) 

As a result, the tight binding Hamiltonian (1) can be diagonalized: 

Ho= Yl e,(k) 7 L(k)7-(k), (6) 

k,<T,!/ = ± 

where the band energies are given by 



e±(k) = e+(k) ± ^6 2 (k) + e%{k) - M . (7) 

In Fig. 2, we show the Fermi surface of this two-band model. The vectors, (0, 7r) and (-7T, 0), 
are the nesting vectors which connect the hole and electron Fermi surfaces. Throughout this 
paper we use the parameters of Rcf. [22]: t\ = —1.0, £2 = 1-3, £3 = t± = —0.85, and \x = 1.45. 
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Fig. 2. The Fermi surface of the two-band model proposed by Raghu et al. The solid (dashed) lines 
represent the hole (electron) Fermi surfaces, and the arrows indicate the nesting vectors. 



By adding the BCS-type interaction Hamiltonian %' , we obtain the following two-band 
BCS Hamiltonian Tl: 

H = Uo + H', (8) 

U> = -J2 [A,7^(k)7i(-k)+A, 7 ^(-k) 7i , t (k)] , (9) 

k,i/=± 

where A u is the superconducting gap of the band v. In this model, the condition A + = — A_ 

(A+ = A_) corresponds to the s_| -wave state (s++-wave state). Using this Hamiltonian (8), 

the normal and anomalous Green's functions in the band v are as follows: 

iu n + £i,(k) 



G u (k, ioj n ) 



Jv )ff ff(k,«*; n ) = Fl aa (k,iuj n ) = - — . , (10) 



where E u (k) = -\/£^(k) + A^ and a =^,4-= 1, — 1. In numerical calculations, we use the 
following interpolation formula: 



A^A^tanh -lj. (11) 

3. Calculation of the response functions 

Now we calculate the response functions for the two-band BCS model. By using the Kubo 
formula, the density response function can be written as 



D(q,iu} m ) 



-J2 ! dTe^ T (T T p r ( q ,r)p t (-q,0)), (12) 

rt J ° 

where r,t denote the orbital indices, and p r (q, r) is defined by 

p r (q,r) = ^4 CT ( k + q> r K CT (k,r). (13) 



kcr 
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Rewriting Eq. (12) in the band representation, we have 

kk' cc' rt vi,-, "4 

x < 1 (k + q)< 2 (k)at 3 (k'-q)at 4 (kO 

x (T r7 t i(j (k + q, r) 7 , 2CT (k, r)^ (k' - q, 0) 7 ^ (k', 0)). (14) 

Now we consider the one-loop contribution to the density response function. Then we 
obtain the following expression: 

D(q,zu; m ) = ^EEEE 



n kcr rt vv' 



(15) 



x a£(k + qK,(k)at,(k)at(k + q) 

x [g u /(k,iuj n )Q u (k + q,iu n - iu m ) 

- -^',a CT (k,iw„)J r ^(k + q,iw„ + iu m )]. 

Here, we substitute the normal and anomalous Green's functions (10) into Eq. (15). Summing 
over the Matsubara frequencies and performing the analytic continuation ioo m — >■ uj + id, we 
finally obtain 

D(a, ^ + *5) = E E E a ^ k + W a "' ( k )°t ( k + q) 



k rt vv 1 
- l-/QMk))-/(i^(k + q)) l 
_ u - E v ,(k) - E u (k + q) + iS 2 

l-/(^(k))-/(^(k + g)) l 
u + E v ,(k)+E„(k + ci) + i6 2 

/(^(k))-/(^(k + q)) 1 



1 + 



1 - 



+ 



u - E u ,{k) + E u (k + q) + iS 2 
/(^(k))-/(^(k + q)) 1 



1 + 



E v ,{k) 

6/(k) 
E u ,(k) 

6/(k) 
E v ,{k) 



i - 



e,(k + q) 



+ 



1 + 
1 + 



K(k + q)/ £</(k)^(k + q) 
g„(k + q) \ , A^/A^ 
£„(k + q 

Uk + q) 



+ 



^/(k)^(k + q) 

A^A^ 

^/(k)^(k + q) 



1 



A^A^ 



^(k + q) 

w + S I/ /(k)-S I/ (k + q) + i<5 2 \^ ^'(k)A^ ^(k + q)7 E u ,{k)E u {k + q) 
where /(-E 1 ) = l/(e^ E + 1) is the Fermi distribution function. The first and second terms in 
the bracket of Eq. (16) are due to the pair creation and annihilation. On the other hand, the 
third and fourth terms are due to the scattering of thermally excited quasiparticles. The factor 
\{- ■ ■ } is the so-called coherence factor. Four coherence factors are included in Eq. (16). 
Next we calculate the spin response function. The spin response function is given by 

X + -(q,^ m ) = -V [dre^r <T T< S+(q, r)^ (-q, 0)), (17) 

rt J ° 

where r, t also denote the orbital indices. Here, S+(q, r) and 5~(q, r) are defined by 

5+(q, r) = (5-(-q, r)) f = E 4 t ( k + T )M k , T )- ( 18 ) 



The spin response function can be calculated in the same way as the density response function. 
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The final expression for the spin response function is as follows. 

(q, w + is) = - a »( k + ^ ( k ) a t ( k + q) 



2 

k rt vv' 



1 - / (^(k)) - / (E u (k + q)) 1 / + £„,(k) ^ ^ £,(k + q) ^ A„/A 



V 



uj-E uf (k)-E u (k + q) + i5 2{\ E u ,(k)J\ E u (k + q)J E u ,(k)E u (k + q 
/ (^(k)) - / (£7„(k + q)) 1 / ^ ^'(k) \f 1+ &,(k + q) ^ A^/Aj, 



w + ^(k)+^(k + q) + i5 2 iV ^'(k)/V ^(k + q)/ ^(k)^(k + q) 
/(^(k))-/(^(k + q)) 1 | &,,(k) \ / | g„(k + g) \ | A«,A„ 



w - £v(k) + E u (k + q) + i<5 2 \ \ E u ,{k) J V E u (k + q) J E u ,(k)E v (k + q) 
+ /(^(k))-/(^(k + g)) If/ Mk)_ \ A &(k + g) \ | A V >A V 



Lu + E„,(k)-E v (k + (i) + i5 2\\ E„,(k) J \ E v (k + q)J E u/ (k)E u (k + q) J _ ' 
Now we present an explanation for the coherence factors. The following two regions are 
considered here: (i) u » |Al 0) | + |aJ?| and T«T C . (ii) w<|Al 0) | + |aJ } | and 0<T<T c . In 
the region (i), the dominant terms in the bracket of each response function are the first and 
second ones (i.e., pair creation and annihilation terms). This is because the third and fourth 
terms vanish due to the fact that / (E u /(k)) — f (E v {k + q)) « for T<CT C . Hence we focus 
on the following coherence factor: 

I Mk) Yl ^ (k + q) V ^ 1 (20) 

2\V ^(k)A ^(k + qV ^(k)^(k + q)J- 1 ' 

The upper (lower) sign in (20) corresponds to the case of the density (spin) response function. 

In the s.) -wave state (A u > A v < 0), the sign in front of the term A v i A v / 'E v i E v in the coherence 

factor at q = Q (nesting vector) is effectively reversed. As already pointed out by Korshunov 

et a/. 23 ) and Maier et al, 24 > 2 ^ this sign reversal causes the resonance peak in inelastic neutron 

scattering measurements. 26 ) In contrast, the same does not hold for the s++-wave state. The 

experiment by Christianson et al. 27 ^ confirms the existence of such a resonance peak and is 

the positive information about the s_| -wave state (see also 5. Discussion). 

Next, we consider the region (ii) on which we particularly focus in this paper. In this region, 

the main terms in the bracket of each response function are the third and fourth ones (i.e., 

quasiparticle scattering terms) because the delta functions 5 (uj — E v i (k) — E v (k + q)) and 

5 (u) + E u /(k) + E u (k + q)) in the first and second terms are zero for almost all temperatures. 

Therefore, we pay attention to the following coherence factor: 

lf(l+M^Yl+ ^ (k + q) ^ 1 (21) 

2 \V i^(k)A ^(k + q)y + ^(k)^(k + q)/- 1 ' 

The upper (lower) sign in (21) corresponds to the case of the density (spin) response function. 

Note that in both the density and spin responses, the sign in front of the term A v i A v /E v iE u 

in the coherence factor in the region (ii) is opposite to that in the region (i). In the region 

(ii), such a sign of the coherence factor directly affects the temperature dependence of the 

response functions. For example, we consider the case A u iA u > (the s ++ -wave state). 
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Fig. 3. Temperature dependence of the density response function normalized by its value at T — T c 
with |A^| = |A^| = 0.05, T c = 0.03, w = 0.001, Q = (tt,0), and 8 = 0.0001. The solid (dashed) 
line represents the case of the s_| -wave state (s ++ -wave state). 




Fig. 4. Temperature dependence of the spin response function normalized by its value at T = T c 
with |A^| = |A^| = 0.05, T c = 0.03, u = 0.001, Q = (tt,0), and 8 = 0.0001. The solid (dashed) 
line represents the case of the s_| -wave state (s ++ -wave state). 



In the density response, the coherence factor cancels the divergence in the DOS, and the 
corresponding response function (i.e., the density response function) decreases monotonically 
as temperature decreases. In contrast, in the spin response, the coherence factor does not 
cancel the divergence in the DOS, and the corresponding response function (i.e., the spin 
response function) has a peak just below T c . In this sense, the coherence factor behaves like a 
"switch" , which controls the contribution from the divergence in the DOS. This is the essential 
point of the coherence effect. 

In Fig. 3, we show the temperature dependence of the density response function. In the 

s_l -wave state, the density response has a peak just below T c . In contrast, in the s++-wave 

state, the density response decreases monotonically as temperature decreases. In the - 

wave state, the gap of the band v and that of the band v 1 have the opposite sign because 
the nesting vector Q(=(7r, 0)) connects the hole and electron Fermi surfaces. Therefore, the 
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Tabic I. The dominant coherence factors in the region (ii) depending on both the time reversal 
symmetry of the external field and the pairing state. We show the case of q = Q = (n, 0). 



even parity odd parity 

(density response function) (spin response function) 



« , wave state i JM X St^i II 1 SijiiH^Li | 4. i"^""" j Q( k ) ) ( 1 4. Mk+Q) \ _ A ,vll A ^-l 1 

s+- wave buaue s ll 1 T b ,ih II 1 ' R..fk+oi I ' k , rki R..fk+oi f 2 ll 1 ' e ,(k) I \ E„(k+Q) J E /(kJE^k+Q)! 



-wave state 




i(k+Q) \ , \A„,\\A V \ 



E u'( k )J\ E„(k+Q) I 1 E^kJE^k+Q)! 



product of the gaps Aj,A„/ has a negative sign, and the upper sign of the coherence factor (21) 
is effectively reversed. As a result, the coherence factor does not cancel the divergence in the 
DOS, and the peak appears. On the other hand, in the s ++ -wave state, the gap of the band v 
and that of the band v' have the same sign. Therefore, the product of the gaps has a positive 
sign, and the upper sign of the coherence factor (21) does not change. As a consequence, the 
coherence factor cancels the divergence in the DOS, and the peak does not appear. 

In Fig. 4, we show the temperature dependence of the spin response function. In the 
s + _-wave state, the spin response shows a monotonic decrease as temperature decreases. In 

contrast, in the s ++ -wave state, the spin response has a peak just below T c . In the s H -wave 

state, the lower sign of the coherence factor (21) is effectively reversed, and the peak does not 
appear. On the other hand, in the s ++ -wave state, the lower sign of the coherence factor (21) 
does not change, and the peak appears. 

Here, we give a brief summary of the coherence effect in the two-band superconductor. 
As mentioned above, it is the sign of the coherence factor that plays a crucial role in the 
coherence effect. In the single-band s-wave superconductor, the sign of the coherence factor 
is determined only by the time reversal symmetry of the external field. However, in the two- 
band superconductor, the effective sign of the coherence factor is determined by both the time 
reversal symmetry of the external field and the symmetry of the superconducting gap (See 
Table I). The coherence factor with the effective plus sign gives the coherence peak. On the 
other hand, the coherence factor with the effective minus sign does not give the coherence 
peak. It is worthy to note that in both pairing states the temperature dependence of the 
density response function is opposite from that of the spin response function. As discussed in 
the next section, this point is quite important for the determination of the pairing symmetry. 

For simplicity, we used the two-band BCS model in our calculation. However, this model 
cannot reproduce the detailed band structure of iron pnictides. This simple two-band model 
and the more realistic five-band model 3 ) can lead to qualitatively different results. For example, 

in the analysis of the superconducting order parameters, Nomura demonstrated that the s H - 

wave state is favored for the five-band model, whereas the p-wave state is for the two-band 
model. 5 ) However, as we mentioned above, the behavior of the density and spin response 
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Fig. 5. Temperature dependence of the density structure factor normalized by its value at T — T c 
with (A^l = |A^| = 0.05, T c = 0.03, w = 0.001, Q = (tt,0), and 8 = 0.0001. The solid (dashed) 
line represents the case of the s_| -wave state (s ++ -wave state). 
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Fig. 6. Temperature dependence of the spin structure factor normalized by its value at T — T c with 
|A^| = |A^| = 0.05, T c = 0.03, uj = 0.001, Q = (vr,0), and 8 = 0.0001. The solid (dashed) line 
represents the case of the -wave state (s ++ -wave state). 





functions in the two-band s-wave superconductor comes from both the peak in the DOS and 
the coherence factor, and it can be simply understood by an extension of the BCS theory for 
the single-band s-wave superconductor. As long as the superconducting order parameter in 
each band has a fully-gapped s-wave symmetry, the details of the band structure does not 
affect the qualitative behavior of both the response functions. Therefore, the qualitatively 
same results are also expected in the five-band BCS model. 

4. Suggested experiment 

In this section, we would like to present an experimental method which determines the 

correct symmetry of the order parameter; the s_| - or wave state. To accomplish this, we 

have to connect the theoretical calculations to the experimental observables. The scattering 
cross section associated with the density (spin) response is proportional to the density (spin) 
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structure factor. The density structure factor <Si)(q, w ) an d the spin structure factor <S x (q, u) 
are given by 

5 D (q,w) = — -^ImDfau + iS), (22) 

7r i — e ^ 

S x (cl,uj) = ~ T -^ J lm X + - (q, u + i8). (23) 

Here, D(q, co+iS) and x H (q,u}+iS) are the density and spin response functions, respectively. 
The temperature dependence of the density (spin) structure factor is presented in Fig. 5 (Fig. 
6). The behavior of the density (spin) structure factor is qualitatively the same as that of the 
density (spin) response function. 

To distinguish between the s_| - and s ++ -wave states, it is necessary to measure the 

temperature dependence of both the density and spin structure factors at the nesting vector 

Q. If the following two behaviors are observed simultaneously, it is expected that the - 

wave state is realized: (i) the density structure factor has a peak just below T c . (ii) the spin 
structure factor decreases monotonically as temperature decreases. On the other hand, if the 
following two behaviors are confirmed simultaneously, it is expected that the s ++ -wave state 
is realized: (i) the density structure factor decreases monotonically as temperature decreases, 
(ii) the spin structure factor has a peak just below T c . 

In real experiments, however, there is no guarantee that the above-mentioned opposite 
behavior of the density and spin structure factors is always obtained. This is because we 
cannot neglect the effect of impurities in iron pnictides. Here we consider the influence of 
impurities on the coherence effect. By using the T-matrix approximation, we calculate the 
effect of orbital-diagonal local impurities 9 ) on the DOS of our two-band BCS model. In Fig. 

7(a) [Fig. 7(b)], we present the DOS in the s_| -wave (s++-wave) state. In the -wave 

state, realistic impurity potentials eliminate the peak in the DOS, and makes the bound state 
inside the gap. Meanwhile, in the s ++ -wave state, the peak in the DOS is quite robust against 
impurities, and the in-gap state hardly appears. If impurities eliminate the peak in the DOS, 
the opposite behavior of the density and spin structure factors is not observed, and both the 
structure factors decrease monotonically as temperature decreases. 

As mentioned in the introduction, the acoustic attenuation and nuclear spin relaxation 
rates represent the density and spin responses, respectively. However, the acoustic attenuation 
measurements observe the density response at small wave vector q ~ 0, and the nuclear spin 
relaxation measurements observe the spin response integrated over q in the Brillouin zone. 
Namely, these two measurements are not suitable for measuring the density and spin responses 
at a large wave vector. Therefore, we hope that the density and spin structure factors at the 
nesting vector Q would be observed by the inelastic neutron scattering measurements. 28 ) In 
the near future, if the resolution is improved, the density structure factor at the nesting vector 
Q may be observed by the inelastic X-ray scattering measurements. 
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Fig. 7. The DOS in the (a) s+_- and (b) s ++ -wave states at T = for n mp = 0.01, \ A^> \ = \ A ( _> | = 
0.05, T c = 0.03, and S = 0.0001. 1 = (solid line), J = 0.5 (dashed line), and I = 4.5 (dotted 
line). 1024 x 1024 k- meshes are used in the numerical calculations. 



At the end of this section we discuss the nuclear spin relaxation rate in iron pnictides. 
In iron pnictides, the coherence peak in the nuclear spin relaxation rate has not been ob- 
served. 29 ~ 33 ) In the s++-wave state, the coherence factor in the spin response at all wave 
vectors in the Brillouin zone does not cancel the peak in the DOS. In addition, in the s++- 
wave state, the peak in the DOS is robust against impurities. Therefore, the absence of the 
coherence peak cannot be explained by the s++-wave state and we can expect that the gap has 

the -wave symmetry. However, in order to determine the pairing symmetry more clearly, 

measuring both the density and spin responses at the single wave vector Q (nesting vector) 
is essential. 

5. Discussion 

In Sec. 3, we mentioned the resonance peak in inelastic neutron scattering measurements. 
There we pointed out that within the one-loop approximation or within the RPA used in Refs. 
[23-25], the peak structure in inelastic neutron spectrum is obtained for ui ~ [A^ | + |A^, | in 

the -wave state, not in the s++-wave state. On the other hand, in a recent study, Onari, 

Kontani and Sato showed that a prominent hump structure in inelastic neutron scattering 
spectrum appears at ui ~ (A^J + |A^ax| in the s ++ -wave state by considering the effect 
of the quasiparticle damping. 34 ) Thus, the peak for oj ~ |A^| + [A^| in inelastic neutron 

scattering experiment may not always the direct evidence for the s_| -wave state. A more 

detailed inelastic neutron scattering measurement at to ~ lA^I + |A^| is needed for the 
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6. Conclusion 

For iron pnictides, we have theoretically proposed an experimental method which distin- 
guishes between the s H - and s ++ -wave states. In each of these two states, we have analyzed 

both the density and spin response functions by using the two-band BCS model within the 

one-loop-approximation. As a consequence, we found that in the -wave (s ++ -wave) states 

the density (spin) response function has a peak just below T c , and the spin (density) response 
function decreases monotonically as temperature decreases. Thus, the temperature depen- 
dence of the response functions in the s H -wave state is qualitatively different from that in 

the s ++ -wave state. We discussed that these results can be explained by the coherence effect 
which depends on the relative phase between the two superconducting gaps. Moreover, from 
the analysis of effect of impurities, the additional information for the determination of the 

pairing symmetry has been obtained. Our suggestion is that the - or s ++ -wave state could 

be identified by measuring the temperature dependence of both the two kinds of response 
functions (two kinds of structure factors in actual experiments) at the nesting vector Q in 
neutron scattering measurements. The theoretically expected behaviors of both the responses 
can be categorized into the following three classes: (i) The density response has a peak just 
below T c , and the spin response monotonously decreases as temperature decreases. In this 
case, we can determine that the gap of this material has the s + _-wave symmetry without 
impurities (i.e., with negligible impurities), (ii) The density response exhibits a monotonic 
decrease as temperature decreases, and the spin response shows a peak just below T c . In this 
case, we can determine that the gap of this material has the s ++ -wave symmetry either with 
or without impurities, (iii) Both the density and spin responses show a monotonic decrease 
as temperature decreases. In this case, we can determine that the gap of this material has the 
s-l -wave symmetry with impurities. 
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